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FORECASTING RAINFALL OF A REGION BY USING FUZZY TIME SERIES 
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Abstract. The definition of fuzzy time series was proposed by Song and Chissom .They presented the time 
invariant model and the time variant model for dealing with forecasting problems. In resent year, many 
method have been proposed to deal with forecasting problems using fuzzy time series .In this paper we 
present a new method to deal with forecasting problems based on fuzzy time series .The proposed 
method belongs to the first order and time variant method Last 18 year total monsoon rainfall (in mm) of 
region Ambikapur, Chhattisgarh are used to illustrate the forecasting process of proposed method. 


INTRODUCTION 


Forecasting activity play an important role in our daily life. The classical time series 
method cannot deal with forecasting problem in which values of time series are 
linguistic terms represented by fuzzy sets .It has the advantage of reducing the 
calculation time and simplifying the calculation process. Chen et al. used the difference 
of the enrollment to present a method of forecast the enrollment of Alabama based on 
fuzzy time series. In [4], Huang extended chen's work presented and used simplify 
calculations with addition of heuristic rules to forecast the enrollment. In [3], chen 
presented a forecasting method based on high -order fuzzy time series for forecasting 
the enrollments of university of Alabama.In [2], chen and hwang presented a method 
based on fuzzy time series to forecast the temperature.[1], presented forecasting first 
order and time variant method for forecasting the enrollment of university of Alabama. 


However the forecasting accuracy rates of the existing fuzzy time series method for 
forecasting enrollment are not good enough. In this paper , we present a method 
belonging to first order and time variant methods It can get a higher forecasting 
accuracy rate for forecasting rainfall of a region Ambikapur in chhattisgarh. 


FUZZY TIME SERIES 


In this section we briefly review the concept of fuzzy time series .The main difference of 
fuzzy time series and traditional time series is that the value of fuzzy time series is 
represented by fuzzy sets rather than real values. Let U be the universe of discourse, 
where U={u1,uz,u3....Un}.A fuzzy sets defined in the universe of discourse U can be 
represented as follows: 


A=fa(u1) /u1 +fa(uz) /U2 +. +fa(un) /Un, (1) 
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where fa denotes the membership function of the fuzzy set A, fa:U-[0,1], and fa(ui) 


denotes the degree of membership of ui belonging to the fuzzy sets A and fa(ui) e 
[0,1],and1<i<n 


Definition :Let Y(t) (t =........0,1,2,......) be the universe of discourse and be a subset of R. 
Assume fi (t) (i =1,2.......) are defined on Y(t), and assume that F(t) is collection of fi(t) , 


f2(t) , f3(t) „u then F(t) is called fuzzy time series definition Y(t) (t=........0,1,2,......). 


Definition: Assume that F(t) is caused by F(t-1)only ,denotes as F(t-1) — F(t) ,then this 
relationship can be expressed F(t) = F(t-1) > F(t) - R(tt-1) , where F(t) = F(t-1) > F(t) 
o R(t,t-1) is called first order model of F(t) , R(t,t-1) is a fuzzy relation between F(t-1) 
and F(t), and © is the Max and Min composition operator. 


Definition: Let R(tt-1) be a first order model of F(t) . If for any t, R(t,t-1) = R(t-1,t-2) , 
F(t) is called time invariant fuzzy time series .Otherwise it is called variant time series. 


Definition :Assume that fuzzified rainfall of i t* year is Ajand that of (i+1) t Ax, where Aj 
and Axare two fuzzy sets defined in the universe of discourse U , then fuzzy logical 
relationship can be represented by Aj—Ax, where Ajis called the current state of the 
logical relationship. 


Method for forecasting rainfall using fuzzy time series 
The proposed method is as follows, it involves following step 


Step 1: Define the universe of U and partition it into several even and equal length 
intervals u1, U2, U3... and Un. For example, assume that the universe of discourse 
U=[500,1900] is partitioned into seven even and equal length intervals u1,u2,u3 u4,Us,U6 
and u7 where 


ui =[500,700] , uz=[700,900], us =[900,1100] , us=[1100,1300] , us =[1300,1500], 
us =[1500,1700] , u7= [1700,1900]. 


Step 2 :Get a statistics of the distribution of the rainfall in each interval. Sort the interval 
based on the number of historical rainfall data in each interval from the highest to the 
lowest .Find the interval having the largest number of rainfall data and divide it into 
four sub -interval of equal length .Find the interval having second largest number of 
rainfall data and divide it into three sub interval of equal length .Find the interval of 
third largest rainfall data and divide it into two sub interval of equal length. Find the 
interval with fourth largest rainfall data and let the length of this interval remain 
unchanged If no data is distributed in an interval , then discard this interval .For 
example , the distributions of rainfall data in different are summarized as shown in 
Table 2. 
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Table 1: Last 18 year total monsoon rain fall(in mm) Data time series region , 


Ambikapur (C.G) 


Rainfall(in 
mm) 


1619.7 

1139.4 
2000 
2002 1086 

1240.6 





Table 2: The distribution of rainfall data. 


Intervals | [500,700] | [700,900] | [900,1100] | [1100,1300] | [1300,1500] | [1500,1700] | [1700,1900] 





Number of 
Rainfall 
data 2 1 5 2 1 1 


After executing this step ,universe of discourse [500,1900] is divided into following 


interval 
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= [500,600] =[600,700] 


u3, = [900, 
u2 = [700,900] 966.66] 
u3,3 = 
ee 66,1033.26] | [1033.26,1100] 


rae 
us2=[1400,1500] 
u7 = [1700,1900] 


step3:Define each fuzzy set Ai based on the re-divided interval and fuzzify the rainfall 





data shown in Table 1, where fuzzy set Ai denote a linguistic value of rainfall data 
represented by a fuzzy setand 1 <i < 14. For example , Ai=very very very very few, A2 
=very very very few ,A3= very very few, A4 = very few, A6s= moderate , A7 = many, As 
= many many , A9 = very many, A10 = too many , A11 = too many many, A12 =too many 
many many, Ai3=too many many many many , Ai4=too many many many many 

many ,defined as follows. 


=1/ u11 + 0.5/ u12 + 0/ u2+0/ u31+ 0/ u324+ 0/ u33 + 0/ u4i+ 0/ u42 +0/ u434+ 0/ 


u44 
+0/ u51 + 0/u52+ 0/u6+0/u7 


A2 =0.5/ u11 + 1/ u12 + 0.5/ u2+0/ u31+ 0/ u32 + 0/ u33+ 0/ u41 + 0/ u42 +0/ u43 + 
0/ u44 


+0/usit+ 0/us2+ 0/u6+0/u7 


A3 =0/ u11 + 0.5/ u12 + 1/ u2 +0.5/ u31+ 0/ u32 + 0/ u33+ 0/ u41 + 0/ u42 +0/ u43 + 
0/ u44 


+0/ u51 + 0/ u52 + 0/u6+0/u7 


A4 =0/ u11 + 0/ u12 + 0.5/ u2 +1/ u31+ 0/ u32 + 0/ u33 + 0/ u41 + 0/ u42 +0/ u43 + 0/ 
u44 


+0/usit+ 0/us2+ 0/u6+0/u7 


As =0/ u11 + 0/ u12 + 0/ u2+0.5/ u31+ 1/ u324+ 0.5/ u33 + 0/ u4i+ 0/ u42 +0/ u43 + 
0/ u44 


+0/ u51 + 0/u52+ 0/u6+0/u7 


Ao =0/ u11 + 0/ u12 + 0/ u2+0/ u31+ 0.5/ u32 + 1/ u33 + 0.5/ u41 + 0/ u42 +0/ u43 + 
0/ u44 
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+0/ u51 + 0/u52+ 0/u6+0/u7 


A7 =0/ u11 + 0/ u12 + 0/ u2+0/ u31+ 0/ u32 + 0.5/ u33 + 1/ u41 + 0.5/ u42 +0/ u43 + 
0/ u44 


+0/ u51 + 0/u52+ 0/u6+0/u7 


As =0/ u11 + 0/ u12 + 0/ u2+0/ u31+ 0/ u32 + 0/ u33 + 0.5/ u41 + 1/ u42 +0.5/ u43 + 
0/ u44 


+0/ u51 + 0/u52+ 0/u6+0/u7 


Ag =0/ u11 + 0/ u12 + 0/ u2 +0/ u31+ 0/ u32+ 0/ u33 + 0/ u41 + 0.5/ u42 +1/ u43 + 
0.5/ U44 


+0/ u51 + 0/u52+ 0/u6+0/u7 


A10 =0/ u11 + 0/ u12 + 0/ u2 +0/ u31+ 0/ u3,2 + 0/ u33 + 0/ u41 + 0/ u42 +0.5/ u43 + 1/ 


u44 
+0.5/ u51 + 0/us2+ 0/ ue+ 0/ u7 


A11 =0/ u11 + 0/ u12 + 0/ u2+0/ u31+ 0/ u32+ 0/ u33 + 0/ u41 + 0/ u42 +0/ u43 + 0.5/ 
u44 


+1/us1 + 0.5/us2+ 0/ u6 + 0/ u7 


A12 =0/ u11 + 0/ u12 + 0/ u2+0/ u31+ 0/ u32 + 0/ u33 + 0/ u41 + 0/ u42 +0/ u43 + 0/ 
u44 


+0.5/ u51 + 1/us2+ 0.5/ uet+ 0/ u7 


A13 =0/ u11 + 0/ u12+ 0/ u2+0/ u31+ 0/ u32+ 0/ u33 + 0/ u41 + 0/ u42 +0/ u43 + 0/ 
u44 


+0/usi+ 0.5/us2+ 1/ Uet+ 0.5/ u7 


A14 =0/ u11 + 0/ u12+ 0/ u2+0/ u31+ 0/ u32+0/ u33 + 0/ u4i+ 0/ u42 +0/ u43 + 0/ 


u44 
+0/ u51 + 0/ u52 + 0.5/ u6 + 1/ u7 


For simplicity, the membership value of fuzzy set Aj either 0,0.5 or 1, where 1 <i<14. 
then, fuzzify the rainfall data in Table 1 and linguistic value of data A1 , Az, A3 psss. A14. 
The reason for fuzzifying rainfall data into fuzzified data is to translate crisp value into 
fuzzy sets get a fuzzy time series. 
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Step 4: Establish fuzzy logical relationships based on the fuzzified rainfall data: 
Aj >Aq, 
Aj Ar 


Where the fuzzy logical relationship" A; A," denotes " if fuzzified rainfall of year n-1 
is Aj ‚then the fuzzified data of year n is Aq ".For example, based on the fuzzify rainfall 
data obtained in step 3, we can get the fuzzy logical relationship as shown in Table 3 


Table3 Fuzzy logical relationship 


ao o 





Step 5: Divide each interval derived in step 2 into four sub interval of equal length, 
where 0.25 points and 0.75 points of each interval is used as upward and down ward 
forecasting points of the forecasting. use the following rules to determine whether the 
trend of forecasting goes up and down and to forecast rainfall. Assume that the fuzzy 
logical relationship is Ai >Aj where Aidenotes fuzzified rainfall of year n-1 and Aj 
denotes fuzzified rainfall of the year n ,then (1) if j > i and the difference of difference of 
rainfall between years n-1 and n-2 and between n-2 and n-3 is positive , then the trend 
of forecasting will go up, and we use Rule 2 to fore cast rainfall ; (2) if j > i and the 
difference of difference of rainfall between years n-1 and n-2 and between n-2 and n-3 
is negative , then the trend of forecasting will go down and we use Rule 3 ; (3) ifj <i and 
the difference of difference of rainfall between years n-1 and n-2 and between n-2 and 
n-3 is positive , then the trend of forecasting will go up, and we use Rule 2 to fore cast 
rainfall ; (4) if j < i and the difference of difference of rainfall between years n-1 and n- 
2 and between n-2 and n-3 is negative , then the trend of forecasting will go down and 
we use Rule 3 to forecast the rainfall ; (5) if j = i and the difference of difference of 
rainfall between years n-1 and n-2 and between n-2 and n-3 is positive , then the trend 
of forecasting will go up, and we use Rule 2 to fore cast rainfall ; (6) if j = i and the 
difference of difference of rainfall between years n-1 and n-2 and between n-2 and n-3 
is negative , then the trend of forecasting will go down and we use Rule 3 to forecast the 
rainfall 


Rule 1, Rule 2, Rule 3 are shown as follows 
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Rule 1: When forecasting the data of year 1997 . there is no data before the rainfall of 
the year 1994, therefore we are not able to calculate difference of rainfall 1995 and 
1994 and the difference of the difference year 1996 and 1995 and between 1995 and 
1994. Therefore ,if |(the difference of rainfall between the year1996 and1995)|/2 >half 
of the length of interval corresponding to the fuzzified rainfall Aj with membership 
value equal tol ,then trend of forecasting of this interval will go upward , and the 
forecasting rainfall falls at the 0.75 of this interval ; if |(the difference of the rainfall 
between the year 1996 and 1995)|/2 = half of the length of interval corresponding to 
the fuzzified rainfall Aj with membership value equal to1, then the forecasting rainfall 
falls at the middle value of this interval ; if |(the difference of the rainfall between the 
year 1996 and 1995)|/2 < half of the length of interval corresponding to the fuzzified 
rainfall Aj with membership value equal to1 ,then trend of forecasting of this interval 
will be downward, and the forecasting rainfall falls at 0.25 point of the interval. 


Rule 2 : If (|difference of differences between years n-1 and n-2 and between n-2 and n- 
3|x 2 +the rainfall of year n-1) or (rainfall of the year n-1 -|the difference of differences 
between yearn-1 and n-2 and n-2 and n-3|x 2) falls in the interval corresponding to the 
fuzzified rainfall Aj with membership value equal to1 ,then trend of forecasting of this 
interval will go upward , and the forecasting rainfall falls at the 0.75 of this interval 
corresponding to the fuzzified rainfall Aj with membership value equal tol, If 
(|difference of differences between years n-1 and n-2 and between n-2 and n-3|/ 2 
+the rainfall of year n-1) or (rainfall of the year n-1 -|the difference of differences 
between year n-1 and n-2 and n-2 and n-3|/ 2) falls in the interval corresponding to the 
fuzzified rainfall Aj with membership value equal to1 ,then trend of forecasting of this 
interval will go downward , and the forecasting rainfall falls at the 0.25 point of this 
interval corresponding to the fuzzified rainfall Aj with membership value equal to1 ; if 
neigther is the case, than we let the forecasting rainfall be the middle value of the 
interval corresponding to the fuzzified rainfall Aj with membership value equal to 1. 


Rule 3 : If (|difference of differences between years n-1 and n-2 and between n-2 and n- 
3|/ 2 +the rainfall of year n-1) or (rainfall of the year n-1 -|the difference of differences 
between year n-1 and n-2 and n-2 and n-3|/ 2) falls in the interval corresponding to the 
fuzzified rainfall Aj with membership value equal to1 ,then trend of forecasting of this 
interval will go downward , and the forecasting rainfall falls at the 0.25 point of this 
interval corresponding to the fuzzified rainfall Aj with membership value equal to1 ; If 
(|difference of differences between years n-1 and n-2 and between n-2 and n-3|x 2 
+the rainfall of year n-1) or (rainfall of the year n-1 -|the difference of differences 
between yearn-1 and n-2 and n-2 and n-3|x 2) falls in the interval corresponding to the 
fuzzified rainfall Aj; with membership value equal to1 ,then trend of forecasting of this 
interval will go upward , and the forecasting rainfall falls at the 0.75 of this interval 
corresponding to the fuzzified rainfall Aj with membership value equal to1 if neither is 
the case, than we let the forecasting rainfall be the middle value of the interval 
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corresponding to the fuzzified rainfall Aj with membership value equal to 1. 
Summarizes the forecasting result of proposed method from 1995 to 2012 


Table 4, Actual rainfall and forecasting rainfall of Ambikapur (C.G) 


Rainfall(in Trend of forecastiig forecasting 
mm) rainfall 


1995 1146.7 
1996 1619.7 middle value 1600 


1997 1139.4 upward, 0.75 points 1137.5 


1998 1049.3 middle value 1066.62 


1225 


1225 
500 
93332 
1066358 
10493 


2008 1358.4 upward, 0.75 points 1375 
2009 603.2 middle value | 6500 


2010 649.7 middle value 
2011 1445.5 middle value 1450 


2012 1181.8 middle value 1175 





In the following, we use the mean square error (MSE) to compare the forecasting 
results with different forecasting method, where the mean square error is calculated as 
follows: 


[ l al 
> (Actual _ra inf all, — forecasted _ ra inf all,) | 
MSE = == 

n 
Where Actual_ rainfalli denotes actual rainfall of year i and forecasting_ rainfall; of 
denotes forecasting rainfall of year i .In following Table 5 we have calculated MSE for 
the proposed method to compare with exiting methods given by [4] and [3] 
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Table 5: MSE calculation 


Year Arthal forecasted rainfall (Ai-Fi)? 
rainfall 


1995 1146.7 
1996 1619.7 1600 388.09 


1997 1139.4 1137.5 
1998 1049.3 1066.62 299.98 
1999 1229.5 1212.5 


1236 1225 121 
2001 1820.5 1800 400 


[2003 s a o 
2004 | e584 | so | 3 
200s | 10668 | 106658 | 00m 
2009 | 6032 | eo | 20024 

oo 


2010 649.7 650 
2011 1445.5 1450 20.25 


2012 1181.8 1175 46.24 
S a EA 





MSE = 2 (Ai - Fi) ° /n =8094.638/ 14 


=578.188 


Conclusion 


In this paper we have investigated using a new method for forecasting the rainfall of the 
region using fuzzy time series. The proposed method belongs to the first order and time 
invariant method.We see that MSE of the forecasting result of proposed method is 
smaller than previous existing methods. That is proposed method gota higher 
forecasting accuracy rate for forecasting rainfall in comparison with the existing 
methods. 
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